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SUMMARY 


A theory is proposed for estimating the noise generated 
at the side edges of part span trailing edge flaps in terms 
of pressure fluctuations measured just in-board of the side 
edge of the upper surface of the flap. Asymptotic formulae 
are developed in the opposite extremes of Lorentz contracted 
acoustic wavelength large/small compared with the chord of 
the flap. Interpolation between these limiting results en- 
ables the field shape and its dependence on subsonic forward 
flight speed to be predicted over the whole frequency range. 
It is shown that the mean width of the side edge gap between 
the flap and the undeflected portion of the airfoil has a 
significant influence on the Intensity of the radiated sound. 
It is estimated that the noise generated at a single side 
edge of a full scale part span flap can exceed that produced 
along the whole of the trailing edge of the flap by 3 dB or 
more . 


1. INTRODUCTION 

A significant proport ion of the noise generated by an 
aircraft on approach to landing is attributable to the deploy- 
ment of high lift airframe components such as slats and flaps. 
Conventional jet-noise is of diminished importance in such 
circumstances because the engines are generally operated at 
reduced thrust. Fink and Schllnker [^], Ahtye, Miller and 
Meecham [5] and Kendall and Ahtye [3] have Identified the 
side edges of part span trailing edge flaps as Important 
sources of the airframe noise. In particular, the extrapola- 
tion to full scale of the experimental results of Kendall and 
Ahtye [5] leads to sound pressure levels which fall only a 
few dB below flight data reported by Healy [^]. Hardin [5] 
has pointed out that the relative intensity of side edge 
sources could be substantially greater for future aircraft 
using high bypass ratio engines, which may well operate in 
the landing configuration with the jet stream impinging 
directly onto the side edges of flaps. 

It is well known (Ffowcs Williams and Hall [^]) that 
unsteady flow in the vicinity of an edge is an important 
source of aerodynamic sound. The generation of noise by tur- 
bulent flow over the trailing edge of an airfoil has been 
extensively discussed in the literature [7]. The studies 
cited above, however, suggest that for deployed part span 
flaps it is predominant ly the characteristics of the unsteady 
flow in the vicinity of the si^de edges which determine the 
level and properties of the edge-generated sound. The unsteady 


motion In the neighborhood of a side edge is related to the 
presence of a component of mean flow, about the edge. This 
flow results in the formation of a conventional tip vortex 
due to forward motion of a lifting airfoil, and is from the 
lower (pressure) to the upper (suction)' surfaces of the flap, 
a separation bubble being formed on the upper surface at the 
edge. The relatively large amplitude unsteady edge flow is 
associated with the instability of this separated flow. 

An idealized two-dimensional theory of the side-edge noise 
mechanism has been discussed by Hardin [5], who considered a 
flow Inhomogenelty in the form of a chordwlse orientated line 
vortex. Potential flow theory was used to model the mean flow 
about the side edge. By means of a detailed study of possible 
trajectories of the vortex, Hardin argued that under a wide 
range of conditions likely to be encountered in practice, the 
noise generated during the motion of the vortex around the 
side edge could be significantly greater than that produced by 
an equal (spanwise) vortex convecting in a mean boundary flow 
past a trailing edge. Such an idealized analysis can only be 
expected to yield possible trends, and, apart from considera- 
tions regarding the characterization of the unsteady edge flow, 
takes no account of Important effects arising from the proxi- 
mity of the wing ahead of the flap and of the side edges of 
the adjacent trailing edge of the wing. The latter are analo- 
gous to the so-called installation effects ^ which play an 
important role in determining the jet-noise of an aircraft 
(see Ref. 8, 9). 

In this report a theory is proposed which permits one to 
estimate the Importance of proximity effects on the side-edge 
noise sources and, in addition, the Influence of varying the 
forward flight speed of the aircraft. To do this the trailing 
edge of the wing is modeled by means of a thin, rigid half- 
plane in the ’’clean" configuration. When a part-span flap 
is deployed the gap between the side-edge of the flap and the 
undeflected portion of the neighboring trailing edge is repre- 
sented by means of a narrow slot in the half-plane which is at 
right angles to the trailing edge. The length of this slot 
is equal to the chord of the flap, and the width is taken to 
be equal to the mean distance between the side edges of the 
flap and the adjacent undeflected wing (see Fig. 1). The 
spanwise edge flow resulting from the additional lift generated 
by the flap is depicted schematically in Fig. 1(c). In using 
this model to calculate the radiated sound, the finite angle 
of deflection of the flap is neglected, except Insofar as 
it determines the magnitude of the mean flow velocity about 
the side edge. This approximation should be adequate for 
moderate angles of deflection, say of the order of 10°, in 
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which case the spanwise component of the mean velocity on 
the upper surface close to the side edge is typically of 
order 0.2-0.3U, U being the forward flight speed [^(?]. 

The analytical problem of calculating the sound generated 
by turbulent fluctuations in the spanwise flow near the side 
edge is formulated in §2 in terms of the evanescent wave 
theory of Chandlramanl [11] and Chase This Involves 

the use of the aeroacoustlc Green’s function whose normal 
derivative vanishes on the slotted half-plane. Approximate 
representations of this Green’s function in the opposite ex- 
tremes in which the chord of the flap is large/small relative 
to the Lorentz contracted acoustic wavelength are obtained in 
Appendices A and B. These are used in §3 to determine the 
spectrum of the acoustic radiation in terms of the wavenumber- 
frequency spectrum of the surface pressure fluctuations just 
in-board of the side edge of the flap. The introduction of 
further simplifying assumptions (§§^, 5 ) enables the acoustic 
spectrum to be expressed directly in terms of a surface pres- 
sure point frequency spectrum. Interpolation between the low 
and high frequency limits in this case yields the approximate 
dependence of the radiation on all relevant Strouhal numbers 
(based on the chord of the flap) and forward flight Mach 
numbers. In §6 a comparison is made with recent experimental 
results^ and it is predicted that the noise generated at the 
side edge of a part span flap at full scale can exceed that 
generated along the whole of the trailing edge of the flap 
by up to 3 dB, 


2. FORMULATION OF THE AERODYNAMIC SOUND PROBLEM IN TERMS 

OF DIFFRACTION THEORY 

A semi-inf Inlt e rigid plate which occupies the half- 
plane (Xj<0, X =0) of a rectangular coordinate system 
(Xi 3 X 25 X 3 ) is taken to model the trailing edge of the "clean" 
airfoil. The aerodynamic sound generated when a part-span 
flap is deployed will be calculated by Introducing a rec- 
tangular slot of length L and width 2s, with 2s << L, to 
represent the side edge gap between the flap and the un- 
deflected trailing edge. The slot occupies the region 
( yL<x j <0 , X 2 = 0 , I X 3 I <s ) of the half-plane, as Illustrated in 
Fig. 1(b). The present analysis will take no account of 
the finite ^ angle of deployment of the flap, although its 
Influence is implicitly Included in that it determines the 
magnitude of. the mean flow about the side edge and the 
associated turbulence levels. In addition, the possible 
existence of a spanwise slot ahead of the flap [indicated 
by the dashed lines in Pigs. 1(a), 1(b)] is Ignored, as is 
the small finite extension of the flap beyond the nominal 
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PART SPAN FLAP 
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UNDEFLECTED 
TRAILING EDGE 


SPANWISE SIDE-EDGE 
FLOW 


FIG. 1. (a) THE "CLEAN" TRAILING EDGE IS REPRESENTED ANALYTICALLY BY 

MEANS OF A SEMI-INFINITE RIGID PLATE IN A MEAN FLOW AT SPEED U 
(b) THE SIDE-EDGE GAP FORMED WHEN A PART-SPAN FLAP OF CHORD L 
IS DEPLOYED IS MODELED BY A RECTANGULAR SLOT WHOSE WIDTH 2s IS 
EQUAL TO THE MEAN WIDTH OF A SIDE-EDGE GAP. (c) VIEW OF THE 
SLOTTED HALF-PLANE FROM DOWNSTREAM. IN CALCULATING THE NOISE 
THE FINITE ANGLE OF DEPLOYMENT OF THE FLAP IS NEGLECTED. THE 
FLAP OCCUPIES THE REGION (-L<Xi<0 ,X2=0 ,X 3 >s) . 
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trailing edge (xi=X2=0) of the wing. Forward flight is 
represented by a uniform mean flow on both sides of the 
airfoil at subsonic speed U in the positive direction of 
the Xj-axls. 

The instability of the separated spanwise flow over the 
upper surface of the flap near the side edge gives rise to 
turbulent surface pressure fluctuations whose magnitudes in 
the immediate vicinity of the edge are influenced by the pre- 
sence of the slot. It is convenient to characterize these 
surface pressures in terms of an hypothetical surface pres- 
sure P(Xi 5X3,1) say, that the same turbulence would Induce 
on a plane rigid wall lying in x^ = 0 . If it could be 
assumed that the turbulence convects in a frozen pattern 
over the flap, P would be equal to the surface pressure 
fluctuation at distances in-board from the edge which exceed 
the correlation scale of the turbulence. The formulation of 
the noise problem in terms of P is equivalent to the use of 
"evanescent wave theory" in the calculation of trailing edge 
noise . The principal advantage in the use of P is 

that it may be assumed to possess local spatial statlonarlty 
in the spanwise direction, a property not shared by the actual 
surface pressure close to the side edge. In practice, surface 
pressure measurements in-board of the edge would have to be 
used to quantify turbulence decay, etc., and thereby to esti- 
mate an approximate form of P near the edge (c.f., the analog- 
ous measurements performed by Brooks and Hodgson . 

Consider an elementary turbulent boundary layer dis- 
turbance above the flap (in X2>0) which is proportional to 
exp {i (k 1 x 1 +k 3 X 3 -ojt ) } . Henceforth, we shall omit from our 
formulae explicit reference to the exponential time factor 
e i^t. Practically all of the turbulence energy will be 
associated with disturbances having subsonic phase velo- 
cities in the ( x 1 -x 3 ) -plane . Hence, if the flap and airfoil 
were absent, the corresponding perturbation pressure would 
decay exponentially with perpendicular distance below the 
boundary layer, and would form an evanescent wave whose 
pressure signature pj could be expressed as 

l(k X +k X )+Px 

p^ = 1 3kg ,m)e (1 ) 

where f = P(k ^ ,kg ,w) > 0 . The functional form of P will be 
discussed below. 
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Just above the flap, at distances from the side and 
trailing edges which exceed ~l/(k )^, the harmonic dis- 
turbance ( 1 ) gives rise to a specularly reflected pressure 
field such that the net surface pressure P(x^,X 25 t) Is 
given by 


= P(k^ ,k^ ,m)e 


l(kjXi+k3X3) 


where , 


P 


= 2 p 


I ' 


( 2 ) 


(3) 


Sound waves are generated when the evanescent wave (1) Is 
diffracted by Inhomogenelt les such as the edges of the flap. 
The side-edge flap noise corresponds to the diffraction radia- 
tion produced at the slot, and use of Eq. (3) In the diffrac- 
tion calculation enables the acoustic pressure to be related 
directly to the surface pressure on the flap In-board of the 
side edge. 

The diffraction problem Is formulated as follows. Let 
pg denote the diffracted component of the perturbation pres- 
sure. The net pressure Is accordingly given by 


P = Pi + P 3 


( 4 ) 


and Is such that on the upper surface of the flap 
Ps Pi , P P 


at large distances from the slot and trailing edge. At large 
distances from the airfoil, Pg represents the acoustic radia- 
tion. In addition, Ps must be chosen to ensure that the 
normal velocity vanishes on the surface of the airfoil. In 
the high Reynolds number limit which Is relevant In practice, 
unsteady viscous stresses may be neglected, and the X 2 -component 
of the momentum equation Implies that the normal component of 
velocity will vanish on the airfoil provided that 
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0 


(5) 


on the upper and lower surfaces of the slotted half-plane. 

Acoustic disturbances are convected by mean flow. Since, 
however, the nonuniform mean flow In the vicinity of the slot 
Is localized. Its effect on propagation will be neglected, and 
In the linearized approximation It therefore follows that p 
satisfies the tlme-harmonlc convected wave equation: 

{(-^*^0 + = 0 . ( 6 ) 


where M = U/c < 1 Is the mean flow (flight) Mach number, 
kp = tjj/c the acoustic wavenumber, and c Is the speed of sound. 
The solution which satisfies the radiation condition of out- 
going waves at large distances together with the surface 
condition (5) can be written down with the aid of a Green's 
function G(x,y). This satisfies dG/dx^ = 0 on the slotted 
half-plane , “the radiation condition, and equation (6) when 
the right hand side is replaced by the point source 6(x-y). 
According to the reverse-flow reciprocal theorem [^5], Grx,y) 
Is also a solution as a function of y of the adjoint or — — 
reverse flow problem: 

{(-i^ - I^) - = 6(x-y) , (7) 

and satisfies the radiation condition and 3G/3y2 = 0 on the 
slotted half-plane. In the reciprocal problem the source 
6(2f-y) Is situated at the observer position ^ and generates 
a disturbance which propagates In y-space in~a medium moving 
at speed U In the negative y ^ -direct Ion . 

A formal representation of the radiated sound is obtained 
by means of Klrchhoff's Integration procedure ([:Z5], Chap. 8). 
In equation (6 ) replace the dummy variable x by y, multiply 
by G(x,y) and Integrate with respect to y over the region of 
space bounded by the airfoil and a large sphere of radius R*, 
say. Similarly, multiply Eq. (7) by Ps(y) and Integrate over 
the same region. Subtract the two results and apply the 
divergence theorem to obtain, as R* «>: 
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(y)dyjdy^ 


( 8 ) 


(x,to) 
s — 


0 G(x,y) 


j 



SP, 


9y 


n 


where the integration Is over the upper and lower surfaces 
respectively of the slotted half-plane, and y^^ Is In the 

direction of the outward normal. In obtaining this expres- 
sion the radiation condition has been used to eliminate 
contributions from the surface of the sphere, and use has 
been made of 3G/3y2=0 on the rigid surface of the airfoil. 
Equation (8) is expressed in terms of the incident^ evanes- 
cent wave Pj by means of equation (5), expressed in the form 

^Ps^^^n "" ^ ^Pl/^y2 ^9) 


on S_^. If we next define 


= G(x,y^ ,+0,y3) 


(10) 


equation (8) becomes: 





[G_j_(x)-G_(x) ] 



(y)dy3dy3 , 


( 11 ) 


the integration being over that portion of the upper surface 
of the slotted airfoil wherein Spq/By^ is different from zero. 
Using equations (1), (3) and summing over all possible values 
of the wavenumbers k^, k^, the radiation is expressed, 
finally, in terms of the surface pressure amplitude coeffi- 
cient ?(k ,k ,03) : 

1 ’ 3 ^ 


p^(x,m) = 


dk^dk 3 T(kj ,kg ,m)P(kj ,k 3 ,(jj) [G^(x)-G_(x) ]< 


i(k^y^+k3y3) 


dy ^dy 

. ( 12 ) 
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Equation (12) may be used to express the frequency 
spectrum S(XaW), say, of the acoustic radiation in terms 
of the wavenumber-frequency spectrum IlCkjjkg^w) of the 
surface pressure PCx^^Xg^t). We recall that these quanti- 
ties are defined through the relations: 


<Pg(x,t )Pg(x,t+x)> 


O/ \ -lOJT, 

S(Xj03)e doj 

J 

^00 


(13) 


<P(yj sy3.t)P(y^+y^,y3+y3 ,t+T)> 

r/T“ 

= n(k,,k .o))e dk dk dw , (1^) 

J J J 1 3 13 

— *oo 


in which the angle brackets denote an ensemble average. Intro 
ducing the hypothesis that P(yi,y 3 ,t) is locally a stationary 
random function of its arguments, with zero mean (a reason- 
able assumption inasmuch as P is defined Independently of the 
diffracted field) , it follows that : 

<P(k^,k3,oj)P(k^,k^,w)> 

= 6(k^+k J6(k3+k3)6(a)+^)n(k^,k^,u3) , (15) 

and similarly that : 

<Pg (x ,ca)Pg (x,oo) > = 6 (a3+a))S(x,oo) . (l6) 

These results will be used below. 


3. THE RADIATED SOUND 

The practical Implementation of the formal solution (12) 
requires a knowledge of the Green’s function G(xjy)* Approxi- 
mate analytical representations of G are available in the two 
asymptotic extremes in which the chord L of the flap satisfies 

k^L/(l-M^) >> , << 1 . (17) 
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The details of the calculation of G In these two cases are 
outlined In the Appendices under the additional restriction 
that kQS'<< 1, l.e., that the acoustic wavelength Is large 
compared with the side edge gap. 

Attention will be confined to the side-edge noise problem 
In which may be taken to be nonzero only In the 

Immediate vicinity of the slot on the upper surface of the 
flap ( -L<y 1 <0 ,y 2 =+ 0 5 y g >s ) . The nonoompaot slot (kpL>>l) Is 
discussed In Appendix A, and It follows from Eq. (A20) that, 
when the observer Is In the far field, and yg>s (on the 
upper surface of the flap). 


G_j_(x)-G_(x) 



( 18 ) 


where 


F 


NC 


-1 

TT^ (l+McosB ) {1 + ^‘InCCk^ s -sln©/^ (l+McosO) ] } 


(19) 


and ^ = 1 . 78107 .., Is the Euler-Mascheronl constant. In these 
expressions the observer position relative to the airfoil Is 
specified In a form appropriate to aircraft flight at speed U 
in the negative x^-direction by spherical polar coordinates 
( r , 0 , (}) ) at the t'lme of emi-ssi-on of the r>eoei-ved sounds where 
the polar angle 6 Is measured from the positive direction of 
the X -axis [see Eqs. (A.I 8 ) of Appendix A]. The emission- 
time (distance r Is given by 


(1~M^ 


-M(x^-y^) 

(1-M^)"'^ 




(x.-yp 


+ xf + X' 


1-M' 


( 20 ) 
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Substituting from (l8) into Eq. (12), we find, as k r^oo 


p^(x,w) 


dk , dk 




13 2r 


0 ■' 


dy^ 


In 


-L 


- 1 

® T 


<iy3 

( 21 ) 


in which the spanwlse Integration has been extended to y ^=+°° 
since in practice the span of the flap will be large compared 
with the characteristic turbulence scale. Convergence of the 
Integral in this limit is ensured by assigning to kg a small 
positive imaginary component which is subsequently allowed to 
vanish. The distance r^ is defined by Eq . (20) when y^ = 0, 

and to leading order we have: 


r = - y ^cos0/(l+M) (1+MCOS0 ) 


( 22 ) 


where use 
from this 


has been made of 
and the integral 


the relations 
definition of 


(Al8). It 
the Hankel 


(z) ([!?"], page 955), that 


follows 

function 


/s (' I 

r p "^(kgS) sin{^L[k ^ -kp cos0/ ( 1+M) ( 1+MCOS0 ) ] } 

p^(x,oj) jj ^^3^0 [k ^ -k|jCos0/( 1+M) (1+MCOS0 ) ] 


X exp{i[k^r^ ~^L(k ^-k^cos0/(l+M) (1+MCOS0 ) ) ] Idk^dkg . (23) 


Let Sj^^(co) denote the frequency spectrum of the acoustic 

pressure in the noncompact limit, defined as in Eq. (16). 
Equation (23) can be used to form the product on the left of 
(l6) which, used in conjunction with Eq. (15), yields: 


11 





I^NcI 


4r' 


rH^^^(k3s) 


k 


nCkj^k^,^) X 


sln^ {^L[k -k cos6/(l+M) (1+Mcos0 ) ]} 

X — — - -— — dk , dk 

{k -k cos0/(l+M) (1+MCOS0) } ^ 

1 0 


( 24 ) 


The formulae given in Appendix B for the Green’s function 
for a aompaot slot [k,j L/( 1-M^ ) <<1 ] may be used similarly to 
calculate the corresponding acoustic pressure spectrum 
say. One finds in this case: 


S(.(co) - 


77 


2 





00 


.oo 


ff (k^L)H 

k 


( 1 ) 

0 


3 


2 

(k3S) I 


n(k 


k3,w)dk^dk3 ,(25) 


in which 


F 


C 


t-T t- 

2k Q L{ (M+cos 0 )^ + (l-M^)cos^c|)sln^0}^ ^ 
17 t(1-M2) (1+McosO) 


sin (^0 ) 

7t( l+McosO ) 


( 26 ) 


f (z) 


1 

m(A)e ^^^dA 

0 


(27) 


The function m(A)(0<A<l) is defined in Appendix B, and depends 
on the value of the slot width/length ratio 2s/L. Its behavior 
for 2s/L = 0.02, 0.05, 0.1 is Illustrated in Fig. 2. The 
angle 0 is given in terms of (0,(j)) by Eq . (B19). When 0 = 0, 
l.e., in the "flyover" plane, we have 

sln(5g0) = sln(5g0){(l-M)/(l+Mcos0)}"^ . (28) 
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m (A.) 


FIG. 



ILLUSTRATING THE VARIATION OF THE FUNCTION m(A) DEFINED IN 
APPENDIX B AND USED IN THE DEFINITION OF THE COMPACT GREEN'S 
FUNCTION, FOR 2s/L = 0.02, 0.05, 0.1. 



It should be noted that Eqs. (24), (25) specify the 
acoustic spectrum in terms of the frequency to as measured 
by an observer fixed relative to the airfoil. This is 
related to the frequency to' that would be measured by an 
observer at rest relative to the mean flow by the usual 
Doppler relation to ' ' = ' to/Cl+Mcos0’) . 


4. APPROXIMATE EVALUATION OF THE WAVENUMBER INTEGRALS 

The integrals in Eqs. (24), (25) can be estimated by 
first noting that, in the region in which the surface pressure 
spectrum H(kj,k 3 ,to) is significantly different from zero, 
k^L >> 1 (I'.e., turbulence scales are small compared with 
the chord' of, the flap),, and therefore that the sine in Eq. 

(24) and f(kiL> in (25) will vary'' rapidly with k^. In the 
case of Eq. (24), this Implies that the principal contribu- 
tion to the kj-lntegral Is'from the neighborhood of 


k* = k^cos0/(l+M) (1+MCOS0 ) 


(29) 


At this relatively small value of k^ the surface pressure 
spectrum is likely to be sharply peaked near kg = mo/Ug, where 
Ug is a velocity which is of the same order as the mean span- 
wise velocity component near the side edge. The value of the 
kg-lntegral in (24) may accordingly be approximated by setting 
kg = oj/Ug in the integrand except in the argument of R. Our 
approximate representation of S (co) is therefore: 

iM U 


Sjjc(“) 





r(k*,to/U^,a))H^^^^(us/U3) ^ 


^00 

J. n(k*,k^ ,m)dk3 . 

-P-CX) 


(30) 


Further reduction of this result is achieved by observing 
that, if it can be assumed that the evanescent wave (1) is 
convected by the mean flow in the same way as acoustic waves, 
we have 

r = {k^+k^-Ck^-Mkj . (31) 
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But M = U /c~0.2M (see Ref. 10) for moderate .subsonic mean 
flow.ftach numbers M. It follows that is generally neg- 
ligible compared with unity, and when this is used in (31) 

with k = k*. it can be seen that in ( 30 ) we may take 
1 1 ^ 

r = |k3| = Icol/Ug . (32) 

It may be remarked that this approximation indicates that the 
possible uncertainty in the precise value of the convection 
velocity to be used in the definition of F is unlikely to be 
significant . 

Similarly, k* may be replaced by zero in the Integrand 
of Eq. (30), since for small spanwise Mach number the 
acoustic wavelength will greatly exceed the characteristic 
turbulence scale. This approximation permits the introduction 
of the point surface spectrum $(m), defined by 


^ ^> ( OJ ) 



— CO 


,(o)dk^ , 


(33) 


where is the characteristic turbulence Integral scale 

in the mean flow direction. This defines the local mean square 
surface pressure P(x^,x^,t) through: 


^oo 

<P^> = J <5(a))do3 . (3^) 


The use of Eq. (32), (33) in Eq . (30) gives our final 
approximation for the acoustic spectrum in the noncompact 
limit : 


^NC 


( 03 ) 





(03S/U3) 


^^(03) . 


(35) 


Introducing the same approximations into the evaluation 
of the wavenumber Integrals in Eq. (25), one obtains in the 
compact limit : 
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J 


(36) 


S(,(io) 


TT^ 1^1 \ ' F I 2 |H^^^(ojs/U ) I 2I(2 s/L)*(m) 


in which 


I = 


m( X) I ^dX . 


(37) 


This integral is a function of the slot width/length ratio 
2s/L, and typical values are given in Table 1: 

TABLE 1 


2s/L 

0.02 

0 . 05 

0.1 

I 

0 . 022 

0.033 

0.046 


5. THE INTERRELATION BETWEEN THE SURFACE PRESSURE SPECTRUM 

AND THE ACOUSTIC SPECTRUM 

Equations (35), (36) provide an approximate relation 
between the acoustic frequency spectrum S(m) and the point 
surface pressure spectrum $(m), although in practice $( 03 ) 
should perhaps be regarded as a chordwise averaged spectrum. 
The chordwise correlation scale is of order U /oo . To 
exhibit an explicit relationship between S(o)) an^ $(( 13 ) we 
shall write 


= aU/o3 , ( 38 ) 

in which a is taken to be constant in the region in which 
^(w) is significantly different from zero. In order of 
magnitude it might be anticipated that a-^U^/U, although its 
precise value will not be required in this section. 

The quantity 

A(o)) - 10 . log^ ^ {S(o))/a$(w) (L/r^ ) ^ } , (39) 
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gives the relation in dB between the sound pressure level at 
frequency w and the surface pressure spectrum. We now proceed 
to discuss the dependence of A(co) on the chordwise Strouhal 
number n = u)L/U for radiation directions which lie in the 
flyover -plane (f) = 0 . 

From Eq. (35) it is seen that the noncompact component 
is given by 




= 10*log^ 


H^^^CfJ(U/U3)(s/L)] 


8ttJ2(1+Mcos6) ^{1+(4/tt^ )ln^ [CfJM(s/L)sine/4(l+Mcos8) ]} 


(40) 

The noncompact approximation requires for its validity that 
coL/c >> 1, ojs/c << 1, so that the outer limits of validity 
of Eq . (40) are delimited by 


1/M < <^ < (L/s)/M 


(^ 1 ) 


The corresponding compact component Ap(fi) Is applicable 

for 


^ << 1/M , 


and use of Eq . 


(28) in the flyover plane gives: 


(42) 


A^(S7) 


10 • log 


1 0 


MslnM%9)I(2s/L) | (U/U , ) ( s/L) ] 

4 (1+M) (1+McosB) ^ 


(43) 


Figures 3? 4, 5 Illustrate the predicted variations in 
A^ (dashed curves) and A^^^ (dotted) with Strouhal number 

for three different radiation directions 0, and for forward 
flight Mach numbers M = 0.1, 0.3j 0.5 respectively. The slot 
width/length ratio 2s/L = 0.1 and U^/U is taken to be 0.2. 
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T 


M = 0.1 
2 s/L = 0.1 


VARIATION OF A(J^) (EQUATION (39)) WITH STROUHAL NUMBER 9 . = ooL/U IN THE FLYOVER 
PLANE (|) = 0, AND FOR M = 0.1, U3/U = 0.2, 2s/L = 0.1: 

.... NONCOMPACT LIMIT OF EQUATION (40); COMPACT LIMIT, EQUATION (39); 

POSSIBLE INTERPOLATION CURVE. THE RADIATION ANGLE 0 IS MEASURED FROM 

THE DOWNSTREAM DIRECTION. 










The solid curves represent possible interpolat ive matchings 
between corresponding asymptotic approximations. There is 
a degree of arbitrariness in the choice of the interpolation 
curves, although it is unlikely that the errors Involved 
will exceed about + 3 dB. The divergence of these curves 
in each figure with decreasing Strouhal number indicates an 
accompanying change in effective acoustic source type. At 
high frequencies the radiation has monopole characteristics, 
modified by forward flight, and is associated with the un- 
steady flux through the slot Induced by the side-edge turbulence. 
At low Strouhal numbers, however, that flux is equivalent to a 
dipole source at the edge of a semi-infinite rigid baffle, and 
the field shape peaks in the forward arc (0 = 180°) with a 
null in the downstream direction. 

The change in source type with Strouhal number can be 
illustrated further by means of the field shape plots of 
Figs. 6, 7, 8, which depict the variation of A(n) with the 
flyover angle 0 for = 0.4, 4, 40, and for the values of 
M, U /U and 2s/L used above. The asymptotic formulae for 
Aq, A]\jq have been used to draw the curves for = 0.4, 40. 

The data points at the Intermediate Strouhal number = 4 
have been obtained by making use of interpolation curves of 
the type shown in Figs. 3-5? and the corresponding field shape 
(dashed curves in Figs. 6-8) must therefore be regarded as 
tentative. The plots exhibit the forward flight amplifica- 
tion of sound radiated ahead of the aircraft (0>9O°), although 
its magnitude is dependent on the Strouhal number Q, It is 
particularly noteworthy that, for nominally equal surface 
pressure fluctuations the semi-baffled dipole radia- 

tion at low Strouhal numbers is typically 10-25 dB greater 
than that from the higher frequency monopole sources. This is 
in contrast to the commonly accepted notions regarding the 
rank ordering of source strength with aeroacoustlc source 
type. The difference arises because the effective monopole/ 
dipole strength is proportional to the Hankel function 
h(1)(o)s/U ), which is small at high Strouhal numbers, but 
logarithmically large as Q 0 . This Influence of the finite 
mean width of the side-edge gap is a proximity effect, similar 
to the so-called installation effects which are known to be 
Important in determining the jet noise of an aircraft [5,5]. 

Its occurrence in the present investigation emphasizes the 
importance of taking account of the neighboring undeflected 
portion of the trailing edge of the airfoil in estimating the 
radiated sound pressure level. 
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THE FIELD SHAPE A(ii) {EQUATION (39)) IN THE FLYOVER PLANE <p = 
FOR Q = 0.4, 4, 40, AND FOR M = 0.1, U3/U = 0.2, 2s/L = 0.1. 
THE CURVE FOR 12 = 4 HAS BEEN OBTAINED BY USE OF INTERPOLATION 
CURVES OF THE TYPE SHOWN IN FIGS. 3-5. 
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6. COMPARISON WITH EXPERIMENT 


The practical application of the formulae derived above 
requires a knowledge of the effective surface pressure 
spectrum $(m). Measured values of $( 00 ) are not yet avail- 
able, although some of its Important characteristics can 
be inferred from the very recent experiments of Miller [20l 
in the 40 x 80 ft^ wind tunnel at the NASA Ames Research 
Center. These involved a l/5~full scale model airfoil of 
6.7 meter semi-span with three flaps of span 3.43 meters 
deployed in tandem. The sound generated at the side edges 
of the flaps at mean flow velocities up to 79 msec ^ was 
found to exceed that from the midspan region by more than 
10 dB, the side edge of the leading flaps being the 
noisiest. The leading flap was deployed at an angle of 
20°, and its geometrical configuration relative to the 
adjacent portion of the airfoil corresponds rather well 
with the idealized model of Fig. 1. We shall here assume 
that the associated values of the surface pressure fluctua- 
tions provide a good approximation for estimating the noise 
generated by the side edges of a single deployed flap. 

At low mean flow Mach numbers, it was found that the 
dominant radiation occurred at sufficiently low frequencies 
that the chord L(~ 20 cms) of the leading flap could be 
regarded as acoustically compact — the field shape 
resembled that of a semi -baffled dipole. Cross-correlation 
measurements of the pressure on the flap Indicated that 
these frequencies satisfied: flap thickness < 2 ttUq/(j^< 2 L, 

where Uq~0.8 U Is the characteristic eddy convection 
velocity. Side edge turbulent velocity fluctuations were 
in the range 10 ~ 25^ of the mean speed U, corresponding 
t o r oot mean square surface pressure fluctuations 

Vpse> say, of order 0.01 — O.O 6 of the mean dynamic pres- 
sure. At the trailing edge of t he flap, the root mean 

square pressure f luctuat ion vp|g , may be identified with 
that beneath a turbulent boundary layer on a flat plate, 
and is roughly O.OO 6 times the mean dynamic pressure [2i], 
These order of magnitude results can be used to compare the 
predicted side edge noise with that generated by turbulent 
flow over the trailing edge of the flap. According to 
Equations ( 78 ), (80) of [7], the acoustic frequency 
spectrum of sound generated by trailing edge tur- 

bulence is given at low Mach numbers, and in the flyover 
plane (({) = 0), by 
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S^^(o}) = • sin^ (?e)<J>(oj) , (44) 

'TT-n 2 


where Lg is the span of the trailing edge wetted by the 
turbulent flow, is the spanwise surface pressure cor- 
relation length and, as before 0 (a)) denotes the surface 
pressure frequency spectrum. This result may be combined 
with the compact approximation formula ( 36 ) for the side 
edge noise in order to estimate the relative Importance 
of the two noise producing mechanisms. To do this, we 
consider a characteristic midrange frequency o) for which 
27rUQ/o)i:5L . At small mean flow Mach numbers the side edge 
generated sound pressure level is predicted to exceed that 
from the trailing edge by an amount ASPL, where 

ASPL = 10 log^jjirML/Lg)(Ji^/<l 3 ) H [ 8 tt ( s/L) (U^/U ^ ) ] 


I( 2 s/L)[p|g/p^g]j (45) 

For a flap of chord L = 20 cms deployed at an angle 
of 20°, the mean value of 2s/L - 0.17. Interpolation from 
Table 1 yields I(2s/L) - 0.053. Taking Uq/U = 0.8, 

U 3 /U = 0 . 3 , then gives 

H^^^^[8it(s/L)(U^/Uj)] ^ ~ 0.112 . 

Finally, we take Z\/ Iz = 6 . 5 . This value is appropriate 
for turbulent boundary layer flow over a flat plate, but 
it is probably of the correct order in the present case 
(c.f. , Substituting into Eq. (45), we obtain: 


ASPL 10 log^j0.38(L/Lg)(p|^/p2^)} . (46) 

Setting L = 20 cms, Lg = 343 cms, we deduce that ASPL lies 

in the range -12 dB to 3.5 dB, when pse/pte varies between 
2 . 8 , 100 in accordance with the experimental values quoted 
above. Thus, the noise generated at a single side edge 
of a part span flap can exceed that generated along the 
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whole of the trailing edge of the flap by more than 3 dB. 
The dimensionless form of Eqs. (45) and (46) Indicates 
that this figure extrapolates to full scale without 
change. 

An alternative expression of this conclusion is 
obtained by taking L = Le in Eq, (46). This gives a 
comparison of the effective acoustic source strengths per 
unit length of side edge/tralllng edge. In this case 
ASPL lies between 0.2 — l6 dB, which is entirely consistent 
with the relative magnitude of 10 dB or more observed 
by Miller [2d] , 


7. CONCLUSIONS 

High and low frequency asymptotic formulae have been 
derived v/hich express the acoustic frequency spectrum in 
terms of the wavenumber-frequency spectrum of surface pres- 
sure fluctuations on the upper surface of a part span flap 
just in-board of the side edge. From a consideration of 
approximate versions of these relations it has been con- 
cluded that : 

1. The radiation has the characteristics of a semi— baffled 
dipole at small values of the Strouhal number o^L/U 
based on the chord of the flap and the forward flight 
speed. At high frequencies the radiation pattern 
resembles that of a monopole source. 

2. The degree of Doppler amplification due to forward 
flight of the aircraft is dependent on the Strouhal 
number wL/U, being characterized respectively by 3, 

2 powers of the Doppler factor 1/(1+Mcos0) as 

ojL/U ^ 0 ,oo. As in the analogous case of noise 
generated at a trailing edge, there exists a weak 
dependence on M of the radiation at 90° to the flight 
path . 

3. The radiation efficiency of side edge noise sources 
is governed by the value of the mean side edge gap 
Strouhal number ms/Ug being larger when ws/Ug is small. 

4. Predictions at intermediate Strouhal numbers can be 
made by interpolation between the low and high fre- 
quency formulae. 
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Recent experiments of Miller [20] Indicate that, at 
least for small values of the flight Mach number, 
the dominant radiation from the side edge resembles 
that of a semi-baffled dipole, the chord of the flap 
being acoustically compact. Theoretical predictions 
extrapolated to full scale suggest that the noise 
generated In this case at a single side edge can 
exceed that generated along the whole of the trailing 
edge of the flap by 3 dB. This figure Is consistent 
with Miller's experimental findings. 



APPENDIX A. THE GREEN'S FUNCTION FOR A NONCOMPACT SLOT 


It is required to solve the reciprocal, time -harmonic 
Eq. (7) ^or k^L >> 1. The width 2s of the slot Is assumed 
to be small compared with the acoustic wavelength. It Is 
sufficient to confine attention to the case In which the 
reciprocal source at x Is In the acoustic far field of the 
slot (|x|>>L), and the field point y Is In the Immediate 
vicinity of the slot. The resulting form of G{x,^) Is 
appropriate for use In Eq. (11), In which the aeroacoustlc 
sources are at points y near the slot and the observer at 
X Is In the far field. 

The definitions 

K = k )** 

0 0 

1 -IK M 

G = -G(l-M^)^e ° 


(X -Y ) 
^ 1 1 


(Al) 


(A2) 


enable Eq . (7) to be reduced to: 


+ 


9Y' 


9y 


9y 


+ KMG = 6(X^-Y J6(x^-y^)6(x3-y^) 


(A3) 


When the slot Is noncompact, such that much of its length 
Is many acoustic wavelengths from either end, the Influence of 
diffraction at the ends can be neglected In calculating the 
behavior of G near the slot. The problem accordingly becomes 
locally two-dimensional, and this _feature is exploited by 
Introducing the Fourier transform of G; 


G 


27T 




— 00 


y 


lk(X.-Y, ) 

,k)e dk . 


(A4) 
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Equation (A3) becomes: 


(tT 7T ^ > (A5) 

\3y^ By" / 

where 

Y = (K^-k^)^ (A 6 ) 

and branch-cuts are taken such that on the real k-axls 
Im(y) > 0 when |k| > I j . 

The particular solution of (A5) due to the reciprocal 
source on the right hand side is 

^ ■ Hj^^Yly-il) , (A7) 

where y = (y ^ 3 '^ » 2. “ (x^^x ). Assume for definiteness that 

X 2 > 0 . When the cylindrical wave (p^ is Incident on the slot, 
and when end-effects are neglected, the diffracted component 
(pjP of (j) in y^ > 0 may be expressed in the form: 

^ ^ (y |y-x* I ) + (yI y I ) , ( |y I >>s) . (A 8 ) 

The first term on the right hand side is the specularly re- 
flected field in the absence of the slot due to an image 
source at x* = (-X 23 X ). When the wavelength of the incident 
wave is large compared with the width of the slot, the un- 
steady flow Induced in the slot gives rise to a monopole 

source of amplitude 3 which is represented by the second term 

on the right . In y <0 the influence of diffraction at the 

trailing edge is Important only at points close to the edge. 
The principal component of the field just below the slot is 
therefore given by 

$ = -gH„^^^Y|y| ) , ( |y|»s) , (A9) 

a monopole equal and opposite in strength to that in Eq . (A8) . 
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At distances from the slo;t which are very much smaller 
than the acoustic wavelength^ cj) can be represented by the 
potential of an essentially incompressible fluid, and we can 
set 




A + B • Re 




(AlO) 


where A, B are constants, and z = y^ + ±y 2 ([15] , §66). The 
values of_3. A, B are determined by matching (AlO) with 
(j) = (j). + (f) given by (A7,8) and with (p given by (A9) in their 
respeit ive^j oint domains of validity. 


__ When the reciprocal source at x is in the far field 
(yl^il >> y is well within a wavelength of the slot, 

we can use the following asymptotic expressions for c|) given 
by (A7,8) and (A9) 


(j) - 2Z + 3{l+(21/7r)ln(CY|y |/2)} (y^ > 0) , (All) 

? - - 3{l+(21/7r)ln(CY|y |/2)} (y^ < 0) , (A12) 

where use has been made of the small argument expansion of 
the Hankel function H(i^([ 2 Z], p. 951), C_ = I. 78107 ... is 
the Euler-Mascheroni constant, and 

Z = -i;(2i/7TY I X I I -I 5 (A13) 

is the asymptotic form of the incident wave . as y|x| ^ 
Equations (All), (A12) are formally valid for^ 

2s << |y| << 1 /|y| 

Note that, although |yI can assume all positive values, it is 
only the interval 0 <|xI<|Kq| which is relevant, wherein the 
incident disturbance constitutes a propagating wave 
(y real) . ^ 

Similarly, for |y| e |z| >> s, the local representation 
(AlO) becomes: 
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(Ali)) 


= A + Bsgn(y^ ) *111(2 I j^|/s) . 

Matching this with expressions (All), (A12) we obtain: 
A = Z 


B = 


2ig 

7T 


-2iZ 

1+ ( 21/7T ) In ( Cys/4 ) 


(A15) 


The substitution of these results into Eq . (AlO), and use of 
(A^) leads to the approximate representation of G(x,y) for 
points y well within an acoustic wavelength of the slot. The 
Integration with respect to k may be evaluated by the method 
of stationary phase when |x| yielding 



in which 


R 


{(Xj-Yj)" 


+ x: 



(A17) 


In the aeroacoustic problem an observer fixed at x would 
be stationary with respect to the airfoil, and would therefore 
be translating at velocity (-U,0,0) relative to the mean flow. 
The practical, flyover, case corresponds to an observer at 
rest relative to the ambient medium, with the aircraft in 
motion at speed U in the negative x ^ -direct ion . It is then 
convenient to express the radiation in terms of emission 
time coordinates (r,0,cf)), say, which give the observer posi- 
tion relative to the source at the time at which the received 
sound was generated. The polar angle 0 is measured from the 
positive direction of the x^-axls, and 
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(Al8) 


Xj-y^ = r(M+cos0) 

X == rcoscbsine 
2 

Xg = rsiri(|)Sln0 
It follows that 

R = r(l+Mcos0)/(l-M2 R - M(Xj-Y^) = r(l-iV[2)^;| 

|x|/R = sln0 ( 1-M^ ) ( 1+Mcos 0 ) , j 

and therefore that expression (Al6) for G(x,y) assumes 
form : 



G(x^y) ^ 


lk„r 
e 0 


4TTr ( 1+Mcos 0 ) 


1 - 



21Re<ln| — + 
s 


I 


- 1 


l+(2i/7T)ln[CkQS*sln6/4( 1+McosG ) ] j 


(A19) 

the 


(A20) 
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APPENDIX B. THE GREEN'S FUNCTION FOR A COMPACT SLOT 

An approximate solution of Eq. (A3) will now be obtained 
in the compact limit for which 

k^L/d-M^) = K„L' << 1 , 

where 

L’ = L/(l-M^)^ . (Bl) 

In this case the length scale of variation of G in the neigh- 
borhood of the slot is small compared with 1/K^ . This implies 
that Kq on the left of Eq. (A3) can be discarded when the 
field point (Y^,y 2 ,y 3 ) is close to the slot, where^ in a lead- 
ing approximation, G will satisfy the Laplace equation. 

When the slot is absent the incident and diffracted 
fields Gq due to the reciprocal source on the right of (A3) 
can be expanded in terms of Bessel and Hankel functions in 
the form: 


Gq = -^ j e^cos[3in(0^-Tr)]Gos[5gn(0-7r) )^n/2^Tp) 


ik(x -y ) 

X e dk (po<p) , (B2) 

where y is given by Eq . (A6), 

Sq = = 2 (n = 1,2,3, . . . ) , 

and (p,0), (pQ, 0g) are polar coordinates defined by 

(X^,x^) = p(cos0,sln0) ; (Y^,y^) = p ^^ ( cos 0 ^ , sln0 ^ ) , 

(see, e.g.. Bowman et at. [i5] Chapter 8). 
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When Kq (X^+x^+x^ ) ^>>la and ( Y^+y^+y ^ ^<<1 , the Integral 
in (B2) can be evaluated by the method of stationary phase, 
leading to: 


^0 


-1 

4irR„ 


00 

I 0 J / 2 (KoPqP/Ro )cos[5gn(F^-Tr) ]cos [^in (^-tt ) ] 
n=0 ' 

l(Kj|R|^-mT/4) 


(B3) 


where 


R„ = (X^+x^+xp*^ . (b4) 

Let Gg denote the additional diffracted field in the 
presence of the slot. This will satisfy Laplace’s equation 
in the vicinity of the slot . At distances from the slot which 
are large compared with its width 2s yet small relative to 
1/Kp (i.e., relative to the acoustic wavelength), Gg can be 
represented by the potential of a line source of strength per 
unit length a(Y^), say, which occupies the portion (-L*<Y^<0) 
of the Y^-axis on the upper surface (y 2 = +0) of the half- 
plane, together with a similarly located line sink of equal 
strength on the lower surface. This can be expressed in 
terms of the incompressible limit of Macdonald’s representa- 
tion of the Green’s function for a rigid half-plane ([iP], 
Chapter 8), from which one finds: 

-sgn(y ) 

^s = J 

-L’ 


o{0 


1 - -- 


tan 


-1 


Rj + (rJ-rO 


d5 . 


(B5) 


where 
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(B6) 


R = {(y^-O" + 

R, = {(Po-O^ + 

The principal branch of tan”^(z)j which vanishes at z = 0, is 
to be taken in (B5). 

Equations (B3), (B5) define 

G = + Gg , (B7) 

In those regions whose closest distance from the slot greatly 
exceeds the slot width 2s, but is much less than the acoustic 
wavelength. Provided that 2s << L, the behavior of G at the 
slot can be taken in a form analogous to (AlO), namely: 

1 z . 

In — + 
s 

where z = ya + ly 23 A Is constant, and B(Yi) varies signi- 
ficantly only over distances which are large compared with s. 

The local source strength B(Y^) Is determined In terms 
of a(Y^) as follows. For fixed Y^ In (-L’,0), and for 
I z I > > s : 

G - A + B( Y^ ) sgn (y 2 ) ln( 2 I z I /s ) . (B9) 

This Is valid In the region 

2s << |z| << L', -L’ < Y^ < 0 . 

The formal expression (B9) must be equivalent to that given by 
(B3, B5j B7 ) as |z| 0 (-L’<Yj<0) In (B7). In this limit 

the Integral in (B5) exhibits a logarithmic singularity of 
the type occurlng in (B9), and equating the coefficients of 
the logarithms we find: 

a(Y^)/TT = B(Y^) . (BIO) 



G = A + B ( Y ^ ) • Re 
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Equations which determine the values of A, B(Yj) are 

obtained by matching the remaining, finite terms in (B7), 

(B9). We do this numerically by partitioning the interval 

-L’ < <0 occupied by the slot into N equal parts, each 

of width h = L’/N. The midpoint of the nth partition is 

at Y, = Y , where 
1 n^ 

Y^ = (n-N-"^)h . (BID 

Within the range (Yj^-^h, Yj^+?gh) B(Yi) is taken to be constant 
and equal to Bj^. A collocation procedure may now be used to 
obtain N+1 equations for A, B^Cn = 1,^, Set Y i = Yj^ 

in Eqs . (B7), (B9). Take the limit |y| e |z| 0 in (B7) and 

match the resulting expression with that on the right of (B9). 
Making use of the identification (BIO) we find: 


A 




47TR. 


(B12) 


B^ln(2h/s) + I hB 


m j IX 

Y -hh " 

m 


" -‘^-5— h - i.tan-' 


) nm 7 t 


Y 

n 


. .1^ 2 


(i?r+iY^i 


0 0 


(2K p Y I /IttR^ ) '^sin (^0 ) (n = 1,2, 


2ttR„ ' 


in which 


(B13) 


e =0 n = m 
nm 


(Bl^ ) 


= 1 n m 


The terms_on the right of (B12), (B13) correspond to the 
inner limit (|y| ^ 0) of which is given by Eq. (B3)* Since 
KjjPq << 1 in the vicinity of the slot, only the first two 
terms in the expansion of in terms of K^p^ have been 
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retained. The component of the Integral on the left of (B13) 
which Involves e^m is evaluated analytically, and that In- 
volving the inverse tangent is approximated by means of the 
trapezoidal- formula. The system of linear equations for 
the Bj^ turns out to be stable provided that the logarithm 
shown explicitly on the left hand side of (B13) exceeds 
unity, l.e. for 2h/s > 2.7. 

Set 


A = |YJ/L' = |yJ/L , (B15) 

and Introduce the dimensionless function m(X) defined in terms 
of B(Yj) by 


B(Yj) 


m(X) 

2ttR„ 


(2K^L»p/17tR^ )^sln(3^e) 


(B16) 


The corresponding discrete values m^^ satisfy Eqs . (B13) with 
the right hand sides replaced by (|Yj^|/L’)^ = { (N-n+%)/N}^, 
(n=l , 2 , • • • ,N) . The local source strength B(Yi) and its 
dimensionless version m(A) are dependent on the slot width/ 
length ratio 2s/L, and the computed variations of m(X) (0<X<1) 
for 2s/L = 0.02, O.O^s 0.1 are depicted in Fig. 2. 

Substituting from (B12), (Bl6) into the local representa- 
tion (B8) of G, using the definition (A2), and expressing the 
final result in terms of the emission time coordinates (Al8), 
we obtain, finally. 


ik r 
0 0 


G(x,y) - 


^Trr „ ( 1+Mcos 6 ) 


+ 2m( y I /L) sin (Jg6 ) [2k L( p/R )/1tt ( 1-M^ ) x 



(B17) 
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in which 


P/Rq = {(M+COS0) ^ + )cos^(j)sln^0} ^/(1+Mcos0) , (Bl8) 


and 

sin (5g0 ) = sgn( )|3s^l-(M+cos0 )/{ (M+cos0 ) ^ + ( 1-M^ ) cos ^ cj)sin^ 0 




(B19) 


Equation (B17) gives the compact form of the Green’s function 
for points y in the immediate vicinity of the slot. 
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